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ABSTRACT 

We present a perturbation theory of kink solutions of discrete Klein-Gordon chains. The un- 
perturbed solutions correspond to the kinks of the adjoint partial differential equation. The 
perturbation theory is based on a reformulation of the discrete chain problem into a partial 
differential equation with spatially modulated mass density. The first order corrections to 
the kink solutions are obtained analytically and are shown to agree with exact numerical 
results. We discuss the problem of calculating the Peierls-Nabarro barrier. 
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I. INTRODUCTION 



In recent years there has been considerable effort in understanding the effects of discrete- 
ness on soliton-like solutions §, §, §, [|, §, §, 0, §, [|, [0]. In this work we will 
restrict ourselves to kink solutions. Kinks connect two groundstates of a choosen system. 
Let us consider a nonlinear Klein-Gordon equation 

To allow for kink solutions the potential V{z) has to have at least two degenerate minima. 
Throughout this paper we will consider only static solutions, i.e. the field $ will be not 
time-dependent. Then equation ( p..l|) is reduced to an ordinary differential equation 

(92$ dV 

^ dx^ ~^ (9$ ^ ' ^^''^^ 
The phase space of ( |1.2|) is two-dimensional. A kink solution corresponds to a heteroclinic 
orbit. This orbit connects the two hyperbolic fixed points (the groundstates) in phase 
space. The invariant manifolds of these fixed points overlap, according to the continuous 
translational symmetry of (|1.1|) , or due to the existence of an integral of motion of (|1.2| ). 



There exist different possibilities to modify the spatial differentials in ( |1.2| ) into differ- 
ences. The most common way is to represent the differences in terms of interaction forces 
between neighbouring particles Xi and Xi^i. 

- + X,_i - 2X1) + ^U=x, = . (1.3) 

Here Xi = $(x = /), and / is an integer (without loss of generality the periodicity of the 
discrete chain is assumed to be equal to one). Equation (|1.3|) is a two-dimensional symplectic 
map, similar to the standard map. In general the invariant manifolds of different fixed points 
do not overlap anymore. Instead they generally intersect in heteroclinic points at nonzero 
angles. The iteration of a heteroclinic point is again a heteroclinic point. One can then 
consider different sequences of heteroclinic points (let us call them heteroclinic orbits). All 
of these orbits will be exponentially attracted to the two fixed points for sufficiently large 



absolute values of /. Exactly two of these orbits correspond to kink solutions, and are thus 
related to their counterparts of the differential equation ( |1.2| ). However these two orbits 
have different energies (in contrast to the differential equation case). The energy difference 
is called Peierls-Nabarro barrier. 

Let us note that there exist also choices of the difference operator such that the invariant 
manifolds still overlap In that nongeneric case static kink solutions exist, which can 
be positioned at any place on the lattice. However, the difference operators are rather 
unphysical, and we will not consider these nongeneric cases here. 

In the limit C — > cxd the two kink-type heteroclinic orbits of ( |1.3|) approach their coun- 
terparts of ( |1.2|) . This is due to the fact, that large values of C imply slow variations of these 
solutions as compared to the lattice spacing. Consequently it is tempting to use a perturba- 
tion approach, which links the kink solutions of (|1.2| ) with the adjoint solutions of ( |1.3|) . In 
this paper we will present a first-order perturbation calculation for the heteroclinic orbits of 
( p..3| ). In section II the difference equation (|1.3| ) is transformed into a differential equation 
with spatially modulated densities. This differential equation is analyzed in section HI with 
the help of separation into slow and fast variables, such that analytical expressions for the 
kink solutions of ( |1.3| ) (in first order perturbation theory) are obtained. In section IV we 
apply our method to two model cases and derive explicit expressions for the kink solutions. 
Section V is devoted to a discussion of the calculation of the Peierls-Nabarro energy. 

II. REFORMULATION OF THE PROBLEM 

Let us consider the following differential equation: 

C<^,,,- p{x)V'{^) = Q . (2.1) 

Here A^^ means (partial) derivative of A with respect to x, and V is the derivative of the 
potential V{z). 

If we choose p{x) = 1, we obtain ( |1.2| ). If we choose 



3 



(2.2) 



l=—oo 



we obtain ( |1.3D |12| ]. This is easy to see by the following reasons. First we note that 
^,xxil < X < (/ + 1)) = or $,^(/ < X < (/ + 1)) = const. Thus it follows 



$(/ + !) -$(/) = <!>,,.(/ + 0.5) 



(2.3) 



By integrating (|2.1| ) from x = I — 0.5 to x = / + 0.5 and using ( p.2|) we obtain 



C{^^S + 0.5) - $,,(/ - 0.5) - V'm)) = . 



(2.4) 



Combining ( p.3|) and (p.4|) we arrive at equation (|1.3| ), where $(/) = Xi. In other words, 
the field $(x) is given by straight lines connecting its values at integer x = I, the field $ 3. 
is given by a function with finite steps at integer x = I and constant elsewhere, and ^^xx is 
a sum over delta-functions, with weights given by ( |2.1| ) using (|2.2| ). 

It is clear that one can make a continuous transition from (3]^) to ( |2.3| ) by varying p{x) 
from p(x) = 1 to (^. 

We rewrite (12. 21) into 



+00 



p{x) = J2 S{x-l) = 1 + 2J2 cos(27rA;x) 

(=— 00 k=l 



Thus we finally arrive at the following equation 



(2.5) 



C^^xx - (1 + 2 ^ cos(27rA;x))\/'($) = . 

k=l 



Note that (|2.6| ) is still an exact reformulation of (|1.3|) 



(2.6) 



III. PERTURBATION APPROACH 



Let us introduce new coordinates x = \fCT and Vt = 27r\/C. Then (12.0) becomes 



$,TT - (1 + 2 ^ cos{kQT))V'{<^) = 
fe=i 



(3.1) 
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In the limit of large values of C the cosine terms in ( |3.1j ) rapidly oscillate due to the increase 
in VL. Thus we can apply standard perturbation treatments using the separation of the field 
$ into slow ^^^^ and fast parts [p!3|: 



$ = $W + ge, . (3.2) 



k=l 



Inserting ( p.2|) into (|3.1| ) and linearizing with respect to the variables S,k we obtain 



$,TT + E ^k,TT = (1 + 2 E COs(mT)) 
A;=l fc=l 



(3.3) 



k=l 

For the fast variables the leading order contribution yields 

efc,TT = 2cos(mT)V'($W) , (3.4) 
ik = -^cos{knT)V'{¥^^) . (3.5) 



Averaging ( p.3|) over the periods of oscillation of the fast variables and using ( |3.5| ) and 
l/e = vrVe it follows 

*!t^t = V:,{^^'^) , (3.6) 

Ve,{z) = V{z)-^{V'{z)f . (3.7) 

Note that equation (|3.6| ) is a simple differential equation, which will be integrated for two 
examples in the following section. 

Since we are interested in the solution of ( p.l|) at integer points, the arguments kQT = 
27m with n being an integer in (|3.5|) . The final solution of (p.l| ) to first order in 1/C is then 
given by 

$(/) = $(^)(/)-^n*^^)) . (3.8) 

Actually ( p.8| ) contains also (incomplete) terms of order One can expand the solution 

in powers of 1/C and extract the first order term after solving along the given path. 
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IV. TWO EXAMPLES 



A. sine Gordon case 



Let us consider 



V{z) 



The kink solution of ( |1.2| ) is given by 



1 — cosf^;) 



(4.1) 



^^"Mx + a) = 4arctan(e^^) 



(4.2) 



Let us consider the slow part of the first-order perturbation. The effective potential 
( p.7|) is given by 

... 1 



sm u 



24C 



sin(2z) . 



(4.3) 



Consequently is the solution of the double sine Gordon equation and can be found in 
jl^ (note that there is an error in equation (3.7) of [|14| - the sign of the power —1/2 has to 
be changed to +1/2) or can be simply calculated by integration: 



^^''\x + a) = 2% - 2arctan 
$^'*'(x + a) = -2arctan 



12C 



cosech 



1 



12Cv^ 



'1 



1 



12C 



-cosech 



'1- 



1 



X 



x > , (4.4) 
X < . (4.5) 



12C 

Here a is an integration constant. Using equation ( pj.8D and expanding in 1/C we finally 
obtain the following first-order perturbation correction for the discrete sine Gordon chain: 

' X + a\ X + a 



$(/) = $^(x + a) + — sech [ — ^ 



2tanh 



(4.6) 



'C J VC 

Since the invariant manifolds of the two relevant fixed points of ( |1.3D do not overlap, but 
only intersect at finite angles, we have to choose the right values of a. Clearly they are 
a = and a = 0.5, which correspond to a kink centered on a lattice site and between two 
lattice sites respectively. These two possible kink solutions are known to exist for the map 

(OD ii- in 



In order to test our result we compute the exact kink solutions of ( |1.3|) with ([4.1|) for 
different values of C. We use the steepest gradient method (minimization of the potential 
energy) and work in quadruple precision. The result will be denoted as Xi. The deviations 
di from its adjoint solution ( |4.2| ) of ( |1.2| ) is then given by di = Xi — $'^(/). The perturbation 
approach yields 0(/) = $(/) — $'^(/) and is defined by the second term on the rhs of ( [4.6[ ). 
In Fig.l we plot di and (pi for C = 10 for both kink solutions (a = and a = 0.5). Clearly 
the perturbation result fits well to the exact one. In order to be more precise, we calculate 
the normalized squared deviation A of the perturbation result from the exact one: 



A 



(4.7) 



2^l=-oo 

Now we can evaluate A for different values of C and see, whether it is monotonously de- 
creasing with increasing C. The results for both kink solutions are shown in Fig.2. No doubt 
the perturbation theory gives the correct first order result. 



B. case 



The second example is given by 



1 



The kink solution of (|1.2| ) is given by 



(^"(x + a) = tanh 



X 



The effective potential ( |3.7| ) reads 



V,s(z) = -(z^-V^ 



2C. 



z 



(4.8) 



(4.9) 



(4.10) 



Thus the slow part is the solution of the differential equation. It can be easily 
integrated using [ |15|] : 



Using (3-8) we finally obtain the first order perturbation result 

<^{^l) = <^(')(^x + a) - -^<l>^'\x + a){<!>^'\x + a)<!>^'\x + a) -1) . (4.12) 

As in the sine Gordon case we calculate di and (pi and plot the results for C = 15 in Fig.3. 
The normalized deviation A[C) is plotted in Fig. 4. Clearly the perturbation theory gives 
the correct result. 



V. THE PEIERLS-NABARRO BARRIER PROBLEM 

Considering the success of the presented perturbation approach with respect to the kink 
solutions, it is tempting to use this result for calculating the Peierls-Nabarro barrier -Epat 
which is given by the energy difference of the two different kink solutions. However as it 
was shown in |jl6[, one has to expect that the leading order asymptotics of EpN contains 



contributions from all orders of the perturbation series for the kink solutions for large values 
of C. This is already clear by noting that the zero-order result (i.e., replacing the exact 
kink solution of the lattice by its counterpart of the adjoint differential equation) yields a 
nonzero Eplf. As shown in [jl6|, these contributions are not enough to fit the exact numbers. 
Clearly at least the first-order perturbation result for the discrete kink has to be taken into 
account (yielding Eplj). But then it is clear that contributions have to expected throughout 
all higher orders of perturbation theory. 

Let us introduce R^^^ = E%/Epn and = E%/Epn, which measure the ratio of 
the first order energy difference (zero order respectively) over the exact one. In Fig. 5 these 
results are plotted for the two examples considered in the previous section. Clearly R^^^ 
is much closer to unity than R^^\ but still there exist discrepancies, which even grow with 
increasing C . This circumstance implies, that the contributions from higher orders of the 
presented perturbation theory in Ep^ gain more weight with increasing C. Consequently 
the first order perturbation scheme, which works as better for the kink solution as larger C 
is, works as better for the Peierls-Nabarro energy as smaller C is. So there is an intermediate 
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range of values of C, where the first order perturbation theory can be satisfactory apphed 
to calculate both kink shape and energy difference. 

Since the calculation of the energy difference using our derived first-order perturbation 
results is still a matter of computing sums, it is not very practically either. Consequently the 
value of the presented approach is clearly in calculating analytically the first-order corrections 
to the kink shape (rather than the energy difference). 

There exists a theoretical approach to obtain the Peierls-Nabarro energy by using prop- 
erties of the splitting angle of the map ( |1.3|) ||T7| . The main point is, that the asymptotic 



dependence of Ep^ on C for large C can be substituted by a form which contains higher or- 
ders. Let us briefly explain this using the sine Gordon chain as an example. The asymptotic 
£'pAr(C)-dependence for large C is given by ||16 



Ep^ = ZCe-'v^(l + 0(i))+0(e^-'v^) . (5.1) 

Here the constant Z is a prefactor, and usually all perturbation approaches are aimed to fix 
this constant. 

To give an example, let us start with the zero-order perturbation with respect to the kink 
shape, i.e., let us take the continuum kink solution ( [1.2D and calculate the energy difference 
Epj^ and then Zc- Naive approaches of this particular task have in addition expanded the 
differences in the energy expression (X/ — X;_i) into series over differentials (see e.g. 0. 
Then noting that the continuum kink shape causes higher order derivatives to decay faster 
with increasing C these higher derivatives were dropped in the mentioned naive approach. 
In the energy expression that means taking into account only squared second derivatives, 
and in the collective coordinate frame (equations of motion for the kink) it would mean 
taking into account only the fourth order derivative of the continuum kink shape |T6[. As 



it was shown in |T^, the dropping of higher order derivatives is wrong, since all of them 



contribute in leading order to Zc- The correct result for Zc gives |16 



2 -| oo 9^m+4^2m+l 

-TT^ - 1 

3 



oo n2m+4^ 

Yi-^r—. -^ 237.82... . (5.2) 

^^2^ ' m(2m + l)! ^ ^ 



To be on the save side we plot in Fig. 6 the ratio 

Z, = E'pj,/{Ce--'^) . (5.3) 

In the hmit 1/C Z^. should approach Z^. Clearly this is the case, demonstrating 
unambiguously that higher order derivatives in the continuum kink shape contribute in 
leading order to Z^- To be precise the naive approach (truncating higher order derivatives) 
yields the result [Q 

647r2 ^ 631.65... 

To compare with the exact result we plot Z for the exact Peierls-Nabarro barrier in Fig. 6 . 
Clearly the zero order perturbation fails by more than 60 percent. 

Finally we plot in Fig. 6 Z for the first-order perturbation result derived in this paper. We 
are about 20 percent off the asymptotic behaviour. Note that both the exact and first order 
result for Z are rather strongly varying with 1/C, so that the true asyptotic dependence 
( p. ID is practically never reached. Here the mentioned approach through the splitting angle 
||T7| has clear advantages. The map (|1.3|) is linearized around one of its hyperbolic fixed 
points. The eigenvalues are given by 



If we replace the exponent (— 7r^\/C*) by (— 7r^/ln(A+)) then we can write 

= ZCe~"'/i°(^+) . (5.5) 

Z is plotted in Fig. 6 for the exact value of Ep^. Clearly the exact result is described by the 
dependence ( ^.5|) very good even down to small values of C, so that: i) ( |5.5|) is describing 
the true Epn{C) dependence much better than ( |5.1|) ; ii) (|0| ) is an easy to use expression; 
iii) the constant Z can be very precisely determined by analyzing the weak dependence of 
Z on 1/C. This yields 

Z ^ 712.26... . 

10 



Having this number, one can determine the exact Peierls-Nabarro energy of the sine Gordon 
chain using (|5.5| ) with very good precision, and down to values C ~ 1. 

VI. CONCLUDING REMARKS 

We have derived first order corrections to the kink shape of a discrete chain. We used the 
methods of slow and fast variables. The resulting differential equations can be integrated 
explicitely, as demonstrated for two examples. Note that the presented method can be 
generalized to the case of anharmonic interactions as well as to time- dependent solutions. 
The generalization of gives 

^5(x-/)$ii-^|,=$,$,,, + ^5(x-/)\/'(<l>) = . (6.1) 

Here W{y) denotes the nearest neighbour interaction on the discrete chain {y = Xi —Xi^i), 
which could well be anharmonic. In the examples considered above we used only harmonic 
interactions 

Wiy) = \cy^ 
so that the second derivative in ( |6.1|) simply yields C. 
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FIGURE CAPTIONS 



Fig.l 

Discrete kink shape deviations from the kink solution of the adjoint differential equation 

versus lattice site I for the sine Gordon chain with C — 10. 

Cirlces - exact result df, 

crosses - first order perturbation result 4>i. 

ei) a — 0.5; b) a = 0. 

Fig.2 

The normalized deviation A of the first order perturbation result (pi from the exact di 
versus C. 

Circles - a — 0.5, crosses - a — 0.5. 
Fig.3 

Same as in Fig.l but for the chain and C — 15. 
Fig.4 

Same as in Fig.2 but for the $^ chain. 
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Fig.5 



The ratio R of the approximated Ep^ over the exact one as a function of C. 

Open circles - zero order result for sine Gordon chain; 

filled circles - first order perturbation result for sine Gordon chain; 

open squares - zero order result for chain; 

filled squares - first order perturbation result for chain. 

Lines are guides to the eye. 



Fig.6 



Z versus 1/C as a test for the asymptotic behaviour of Ep^ for the sine Gordon chain 
(see text). 

Open circles - zero order result; 

filled circles - first order result; 

open squares - exact result; 

filled squares - Z for the exact result (see text) . 

Solid hues are guides to the eye. 

Dashed line - predicted value Z a; 712.26... for the exact result. Clearly Z is approaching 
this value for 1/C ^ 0. Also the exact result for Z (open squares) is approaching the dashed 
line in the same limit. 

Dotted line - Z = 64%'^ 631.65... (see text). 
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